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We study the behavior of exciton polaritons in an optical microcavity with an embedded semicon¬ 
ductor quantum well. We use two-component exciton-photon approach formulated in terms of path 
integral formalism. In order to describe spatial distributions of the exciton and photon condensate 
densities, the two coupled equations of the Gross-Pitaevskii type are derived. For a homogeneous 
system, we find the noncondensate photon and exciton spectra, calculate the coefficients of transfor¬ 
mation from the exciton-photon basis to the lower-upper polariton basis, and obtain the exciton and 
photon occupation numbers of the lower and upper polariton branches for nonzero temperatures. 
For an inhomogeneous system, the set of coupled equations of the Bogoliubov-de-Gennes type is 
derived. The equations govern the spectra and spatial distributions of noncondensate photons and 
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I. INTRODUCTION 

In the last two decades there has been extensive 
work on Bose-Einstein condensation in the system of 
two-dimensional exciton polaritons in optical micro- 
cavites (see the reviewsir— and the references therein). 
The profound investigation of properties of the polari¬ 
ton system has taken place^^— . Bose condensation 
has been demonstratedi^^— and many features of the 
polariton Bose condensate such as superfluidity22,i2^, 
vortices^i^— and solitons^^— formation were studied, 
and have been shown to have specific character in the 
non-equilibrium spatially inhomogenious system of exci¬ 
ton polaritons. Recently, novel phenomena such as bal¬ 
listic transport over macroscopic distance o^^d^ , sponta¬ 
neous oscillations^^— , Josephson-like effects^^— , op¬ 
tical Aharonov-Bohm^ and spin Hall effects^ are ac¬ 
tively discussed. 

Most authors of theoretical works use, explicitly or 
not, the conventional approach consisting of the follow¬ 
ing steps: one writes out the Hamiltonian of the system 
including exciton and photon creation and annihilation 
operators, performs the Hopfield canonical transforma¬ 
tion from the exciton-photon to the lower-upper polari¬ 
ton basis, neglects the upper polariton branch as empty 
at low temperatures, and regards polaritons on the lower 
polariton branch as a dilute Bose gas weakly interacting 
via the exciton component. Such an approach allowing 
the straightforward application of the Bogoliubov the¬ 
ory has led to a great amount of important results (see 
the broad review of the approach in Ref4). 

However, this approach has two serious shortages. 
The first one is the neglection of the upper polariton 
branch. It is quite justified at the temperature range of 
a few degrees Kelvin which takes place in up-to-date ex¬ 
periments. Though, as many authors have noted, Bose 
condensation of polaritons is so attractive due to ex¬ 
tremely high transition temperature being the effect of 
very small polariton effective mass. Already for nitrogen 
temperatures, the upper polariton branch occupation is 


not fully negligible. Second and more serious shortage of 
the conventional approach is that the canonical transfor¬ 
mation from the photon-exciton to the polariton basis 
is performed in the momentum space and hence it re¬ 
quires the spatial uniformity or quasi-uniformity of the 
system. 

With the up-to-date traps having dimensions of the 
order of dozen microns^^— , one can perform the canon¬ 
ical transformation at each point in space in the frame 
of the quasi-classical approximation. However, latest 
experiments have studied systems already with micron 
dimensional^— and it is likely that the trend of minia¬ 
turization will go further. Besides, the Bose condensed 
system can exhibit topological defects such as vortices 
and sohtona2i“— . The exciton and photon healing 
lengths, which determine the spatial scales of these de¬ 
fects, are of the order of micron. Thus, starting from the 
micron scale the conventional approach becomes non- 
applicable. Moreover, polariton as a well defined quasi¬ 
particle at this scale does not make sense. 

In this paper, we study the behavior of exciton po¬ 
laritons in the framework of a two-component exciton- 
photon approach. The exciton and photon components 
are described as independent quantum fields interact¬ 
ing via Rabi splitting. Bose condensation leads to the 
appearance of the two coupled exciton and photon con¬ 
densates. In the mean field approximation this approach 
has been used, for example, in Refs— . The au¬ 
thors of Refs— have gone beyond the mean field ap¬ 
proximation and calculated the energy spectrum of non¬ 
condensate photons and excitons. 

We analyze the equilibrium Bose condensed polari¬ 
ton system without pumping and photon decay from 
the cavity at nonzero, but sufficiently small tempera¬ 
ture. These effects were studied in the two-component 
approach in Ref—. Moreover, we do not take into ac¬ 
count the spin degree of freedom of polaritons. 

The approach is formulated in terms of path inte¬ 
gral formalism for the following reasons. First, this for¬ 
malism allows the study of both condensate and non- 
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condensate parts of the system from the unified point 
of view. The second advantage is the use of coherent 
states, which are the natural choice for an exploration 
of quantum properties of light emitted from the polari- 
ton system. After all, the method involves Matsubara 
technique and hence it is proved to be effective for a 
study of the polariton behavior at nonzero tempera¬ 
tures. The application of the path integral method to a 
weakly interacting condensed Bose gas can be found in 
the excellent review by Stoof^. 

The paper is organized as follows. In Section II we 
describe the geometry of the polariton system and in¬ 
troduce the exciton and photon one-particle eigenstates. 
In Section III we write the partition function of the sys¬ 
tem as a path integral, which contains an action fully 
describing our system. The exciton and photon quan¬ 
tum fields in the path integral formalism are c-number 
fields written in the basis of coherent states. In Section 
IV we turn to the Bose condensed polariton system and 
write out both exciton and photon fields as a sum of the 
order parameter and fluctuations around it. The part 
of the action, which does not contain fluctuations, is 
the Pitaevskii-like functional. Its minimization leads to 
the two coupled equations of the Gross-Pitaevskii type. 
These equations give the description of the system of the 
two coupled condensates in the mean field approxima¬ 
tion. We analyze the properties of the obtained equa¬ 
tions and find the conditions of equivalence between the 
one-component and two-component approaches. In or¬ 
der to explore the noncondensate part of the system, in 
Section V we analyze the quadratic in fluctuations part 
of the action. The most important result here is the 
canonical transformation from the exciton-photon basis 
to the lower-upper polariton basis. This transforma¬ 
tion incorporates both the Hopfield polariton transfor¬ 
mation and the Bogoliubov transformation for weakly- 
interacting Bose gas. Moreover, we obtain nonconden¬ 
sate particles energy spectra and calculate exciton and 
photon occupation numbers for lower and upper polari¬ 
ton branches at nonzero temperatures. In Section VI we 
derive the set of coupled equations of the Bogoliubov-de- 
Gennes type. This equations govern the energy spectra 
and spatial distributions of the noncondensate excitons 
and photons in a spatially inhomogeneous system. 


II. GEOMETRY AND ONE-PARTICLE 
EIGENSTATES 

We study a system of exciton polaritons in a semicon¬ 
ductor optical microcavity with an embedded quantum 
well. In the simple case when the quantum well pos¬ 
sesses in-plane translational symmetry, energy spectrum 
of excitons in the region of small in-plane momenta has 
the form 


do) 


2mp 


( 1 ) 


where e is the dielectric constant of the medium, = 

2TOea:e^/e^?i^ is the 2D exciton binding energy, mex = 
rrie + mh is the 2D exciton mass, and rUe and ruh are the 
effective masses of an electron and a hole, respectively. 

An external confining potential V'^“(x) for excitons 
in the quantum well can be created. Experimentally, 
there are several methods to realize the confinement of 
excitons. One of them is the exciton energy shifting 
using a stress-induced band-gap shif t^^i^^'^^ . 

We assume one-particle eigenstates describing non¬ 
interacting excitons confined by the external poten¬ 
tial V®^(x) to be found from the time-independent 
Schrddinger equation 

The set of solutions Xn(x) is ortho normalized 

J dx X„(x)Xm(x) = (5„m (3) 

and full 

5Zxn(x)x„(x') =(5(x-x'). (4) 

n 

In the case of the planar microcavity the photons in 
the region of small in-plane momenta have the following 
dispersion: 



Here niph = 'kTi^/cL is the effective photon mass. We 
consider the lowest state n = 1. 

There are several experimental approaches to realize 
the confinement of photons^. It is, for example, cre¬ 
ating a special dielectric permittivity profile inside the 
cavity which would lead to photon localization, or mak¬ 
ing a trap for photons based on special shaping of the 
microcavity width T(x) (up to creation of a finite system 
limited by mirrors from all sides). In this case the effec¬ 
tive photon mass ruph and the first term in ([S]) become 
functions of coordinates. As a result, the following one- 
particle problem arises for photons in the microcavity 
with a nonconstant width T(x): 


+ (i(^ - 

We will assume microcavities to possess in-plane trans¬ 
lational symmetry, and will not consider microcavities 
with photon confinement, i.e. TOp/j(x) = niph = const. 
The one-particle problem for photons reduces to 


V'n(x) = 0. 


(7) 
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Similarly to the above exciton one-particle eigenstates, 
we assume the set of solutions of this problem to be 
orthonormalized and full 

J dx Vin (x)'(/'m (x) = (5nm, (8) 

^t^n(x)V!„(x') =d(x-x'). (9) 


III. ACTION AND FREE-FIELD GREEN’S 
FUNCTIONS 


with the coefficients of expansion Xtiir) and ipniT) de¬ 
pendent on time. 

We include the detuning Eq between the exciton spec¬ 
trum = Eo + p'^j2mex and photon one = 
p'^/2mph into Sex and consider iHp > 0. 

Further we will assume the exciton interaction to have 
the approximate form U(x — x') = Vq d(x — x'). 

The partition function Sex [x, x] can be rewritten as 


Sex[x,X] 


r r 

dr dx dr' / dx' 


X(x, t)G (x, t; x', t')x(x', t'), (17) 


In the functional approach, the grand-canonical func¬ 
tion of the exciton-photon system is the functional in¬ 
tegral 

Z = j DxDxDiIj (10) 

where the action S consists of four terms 

s[Zx,'4’,i^] = Sexix.x] + Sph[i^,'4’] 

+ Srabi[x,X.' 4’,'4’] + Sint [x,x], (11) 

describing excitons, photons, Rabi splitting and interex- 
citon interaction, respectively: 


where 

G 


(0)^ 


(xr-x' r') = --^--^ 

’ ’ ’ h \ dr 2mex 


+V^'^{x) + Eq - fi]S{x- x')6{t - t') (18) 


or, equivalently. 


d 71^X7^ 


X Ge,r( 0 )(x,T;x',T') = -M(x-x')(5(r-r'). (19) 


rhP e 

SexlZx] = J dr J dx x(x,r) 

X + W“(x)-f Go-Ai ) x(x,t), (12) 


Sphiip^tp] = J dr J dx ip(x,T) 

eUP e 

SRABi[x,X,i’,tp]= J dr J dx 

X ^ ['0(x, r)x(x, r) -f V’(x, t)x(x, t)] , (14) 


1 f 

S^nt[x,x] = / dr dxdx' x(x,t)x(x',t) 

X U(x- x')x(x',t)x(x,t). (15) 

Here x(x, t) and ^(x, r) are the field operators written 
in the basis of coherent states and hence they constitute 
c-number functions connected with Xn(x) n via 

relations 

X(X, t) = y" Xn{T) Xn(x), '0(x, t) = V V’„(t) V’n(x) 


(16) 


Substituting (5(x — x') from (|3]) and using the equality 
<J(r-T') = ^Ee—(20) 

we obtain 


1 




XXn(x)Xn(x')- — -. (21) 


hp 


The solution of the equation (1191) reads 

-h 


Gex( 0 )(x,T;x',T') = ^ 


-ihujm + efp + Eq - fi 


X Xn(x)Xn(x')- 


^ — iuJrn{T — T') 


np 


( 22 ) 


Since the functions and Gea,(o) ^^re mutually in¬ 


verse, they obey the relation Gea;(o)Gg^^p^ = /, where I 
is the unity matrix. In coordinate space, this relation 
has the form 


dr” J dx" Gea;(o)(x,T;x",T") 


X G-ig)(x",r";x',r') =<5(x-x')<5(r-r'). (23) 
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Similarly, 5^/1 [■0,'0] can be written as 

X 'i/'(x,T)Gp^^(0)(x,r;x',r')V’(x',T'), (24) 


Substituting the definitions (l30ll into the action given 
by m and rearranging the derived terms according 
to their orders in x'(x, r) and tjj'(x,T), we obtain the 
following results. 

The zeroth-order part of the action is given by 


where 


Gp,!(o)(x,T;x',r') = 

“ “ 17: “^ 

This equation is equivalent to 

— hS{x — k')S{t — t'). (26) 

For the functions Gp;i(o) (x, t; x', r') and 
^phio) 

<^17(0) .'r') =-\ - p) 


X ^/;„(x)V’n(x') 




np 


(27) 


Gp/i(o)(x,T;x',r') = ^ 


-h 


-ihiOm + - fj. 


X V'n(x)'i/'n(x') 




np 


(28) 


fhfi p 

J dr" j dx" Gp^(o)(x,r;x",r") 


X G-i (x",r";x',r') =<5(x-x')d(r-r'). (29) 


IV. COUPLED GROSS-PITAEVSKII-LIKE 
EQUATIONS 

We now turn to the Bose condensed exciton-photon 
system. Let us introduce the time independent order 
parameters Xo(x) and tpoi^) S'S 

x(x,t) = xo(x)+x'(x,t), '!/'(x,t) = ?/’o(x)+V’'(x,t), 

(30) 

where x7x, t) and are the fluctuations of the 

exciton and photon quantum fields around the order pa¬ 
rameters Xo(x) and ^/'o(x). We shall use the Bogoliubov 
approximation assuming x7x, t) ^ Xo(x) ■(/’'(x, r) <C 
■ 00 (x), i.e. we confine ourselves to the case of low tem¬ 
peratures. 


?i/37^p[xo,Xo;V'o,7o] = J |^;))~IVXo(x)p 
-b (U®"^(x) +Eo-fi) |xo(x)|^ -b y |xo(x)|'‘ 

+ - mI</>o(x)P 

Zruph 

+ ^ [7o(x)xo(x) -b V'o(x)xo(x)] I (31) 


and has the form of the Pitaevskii functional. 

According to the semiclassical method, the functional 
reaches its minimum at functions Xo(x) and ^o(x), for 
which the functional variation vanishes. Equivalently, 
we can set the part of the action linear in 
^'(x, t) equal to zero. As a result, we obtain the set 
of coupled stationary equations analogous to the Gross- 
Pitaevskii equation: 


2^2 


^ V 

2top 


+ U“(x)+£;o-p 
+ K)|xo(x)P^Xo(x) -b ^V'o(x) = 0, 

7’o(x)-b ^Xo(x) = 0. (32) 


These equations provide a minimum of the functional 
m and describe the system of the two coupled con¬ 
densates in the mean field approximation. 

Let us discuss the derived equations by considering 
first the homogeneous case V“(x) = 0. In the absence 
of the Rabi splitting tiil = 0, the chemical potential 
of excitons equals fi = Eq + Vong®, same as for the 
conventional one-component Bose condensed gas. Here 
we introduced the exciton condensate density |xoP = 
ng®. It is known that the Rabi splitting leads to the 
appearance of the two branches in the energy spectrum, 
which are referred to as the lower and upper polariton 
states (below we denote them by superscripts “L” and 
“U”). The chemical potential is equal to 


d 


iL,U) _ 


1 


(Ao + Ro ng") T ^{Eo + Vml^ f 

(33) 

and varies within the interval G £ 

[Mmil.oo], where ^7(7^ = (^o T \JEl \- h?Q?)/2. 

In both cases the solution of the set of equations (1321) 
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has the form 


Xo 


ipo 


1 /A 


Va 




nn /1 ^2 _ Eon-n^xi^/A 
2/i V vb ^ 


(34) 


For the lower polariton state, the exciton and photon 
condensate phases differ by tt, since the chemical poten¬ 
tial is always negative. For the upper polariton state, 
the exciton and photon condensate phases coincide and 
the chemical potential is positive. 

The energy of the homogeneous system per unit area, 
as it follows from (ED), can be calculated as follows: 


= Eon,, + ^Vo{n^o" f T (35) 


where = |'0oP is the photon condensate density. 

Thus, when the homogeneous exciton-photon system 
is the condensate of lower polaritons, it has the mini¬ 
mal energy E^^\ When the system is the condensate 
of upper polaritons, it has the maximal energy E^^\ 
The difference between these energies per one particle 
is the energy of Rabi splitting fiQ, (in the limit of low 
density). The energies of the uniform states, for which 
phases of the exciton and photon condensates differ by a 
value ranging from 0 to tt, are situated within the inter¬ 
val [E^^\ E'^^'^]. Such states are time-dependent, since 
the exciton and photon condensates cyclically turn into 
each other, i.e. a kind of the internal Josephson effect 
takes place. We shall discuss these states elsewhere^. 
Here we confine ourselves only to the study of the lower 
polariton branch with the energy E^^\ omitting further 
the superscript “L” for shortness. 

Substituting Eq. (IM)l into p5l) . we get the connection 
between the energy and the chemical potential 


E = 


2Vo 


9 1 .. 9^9 9 h^fi^ 

y? A- 2^^^^ ~ -®o- 


0 


We introduce the polariton condensate 
sum 


3 

16 ^2 

(36) 

density as the 


(37) 

Substituting Eq. p4l) into ED, we find the connec¬ 
tion between the polariton condensate density and the 
chemical potential 


n 


p 

0 


1 


1 

fj,- Eq - -——Eo 
4 


1 

16 


(38) 


With the help of Eqs. (15^ and (1551) we can easily prove 
that dE/dno = /i, i.e. the chemical potential is equal to 
the change in the energy as one extra polariton is added 
to the system. The fact that ^ < 0 means that the en¬ 
ergy always decreases. It is worth pointing out that 


the expressions (1341) could be obtained from the condi¬ 
tion that the energy has a minimum at a fixed value of 
chemical potential (clE/Sn^)^ = {dE/dn^)^ = 0. As 
it is seen from (IM)) , in the limit of large density /r —0 
the polariton condensate becomes mainly photonic one, 
and the energy is now connected with the density by 
the relation E = —(3/16)(?i^H"‘/Vo)^/^(n^)^/^, which 
differs principally from the relation for the conventional 
one-component Bose condensate E = (l/2)Vb(u“)^. 

Let us next discuss the inhomogeneous case. The 
lower equation in (1321) . describing the photon conden¬ 
sate, possesses the spatial scale which plays a role of 
the healing length and equals 


^ph — 


\l 2jnpfifl 


(39) 


The length ^ph has a minimal value in the limit of zero 
density and increases with the chemical potential. 

For the upper equation in (1321) , describing the exciton 
condensate, the healing length equals 


^ex 


h 

sj2mex{Eo - m ) 


(40) 


and increases with the chemical potential as well. 

It is important to point out that the relation 
nriph/nn,, ^ 10“^ observed typically in experiments 
causes the following inequality 


^ph ^ C ex ■ 


(41) 


Let’s now study the conditions of the equivalence of 
the two-component exciton-photon and one-component 
polariton approaches. These approaches are obvi¬ 
ously equivalent in the homogeneous case, when '00 = 
{fiQ/2fi)xo, no = |0oP + IxoP- We expect these local 
relations to be valid also in the case of smooth densities 
change. We rewrite the lower equation in (1551) in the 
integral form using the fact that the MacDonald func¬ 
tion Ko{x) (see, for example;^) is the Green’s function 
of the equation. Thus we have 

V’o(i-) = -^^ / (42) 

Here we introduced the dimensionless vector r = x/^ph- 
Taking the asympotics Ko{x) ~ sJ-kI 2xe~^ valid for 
X —>■ oo into account, we can expand the smooth func¬ 
tion x(rO in n power series in the vicinity of r' = r and 
keep only the first terms 

Xo(r') ~ x)!-) + Vxo(r) ' ^ + 

ij * ^ 


where s = r' — r. We substitute this expansion into (1421) 
and integrate it with respect to ds = s ds dip. Here it is 
convenient to use the integral formula for the MacDon¬ 
ald function 
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As a result, the gradient term and the term contain¬ 
ing the second mixed partial derivative vanish after the 
integration with respect to the angle dip, and we obtain 


V-oCx) = ^Xo(x) 


, ^2 Axo(x)' 


(45) 


where A is the Laplace operator. Therefore the condi¬ 
tion that determines the applicability of the local rela¬ 
tions between particles densities is 


Axo(x) J_ 
Xo(x) 


(46) 


which means that the relative exciton density change 
must be small on the scales comparable to the photon 
healing length. 

Ill topological structures such as vortices and solitons, 
the exciton density changes considerably on the scales 
of the exciton healing length. The same picture takes 
place near the boundaries of the condensate confined by 
the box potential. In all these cases the one-component 
polariton approach is not applicable. 

The extensively used Thomas-Fermi approximation 
consists ill neglecting the quantum pressure terms in 
the equations (15^ at sufficiently high densities of con¬ 
densate particles. The obtained estimate establishes the 
range of validity of Thomas-Fermi approximation. 

The solution of Eq. (15^ for the harmonic trap is 
published in Refi^, the vortex solution is discussed in 
Ref,^. 


V. NONCONDENSATE PARTICLES 

A. Matrix Green’s function 

We can represent the quadratic part of the action 
gquad g quadratic form 


gquad = --J^ dr j • G-' 


(47) 


where the Green’s function has a matrix structure 

/-I 


-G = 


X 

V’' 

x' 

L^'J 


• [x'^V^'^x'^V”'] 


(48) 


and can be written as 


G 1(x,t;x',t') = G,„s(x,r;x',r')- 


(0)^ 

'2Ro|xo(x)P Tin/2 yoXo(x)^ 0 

nn/2 0 0 0 

K)Xo(x)^ 0 2Ro|xo(x)P fiQ./2 

0 0 riVi/2 0 

X (5(x — x')(5(t — r^), (49) 


where 


GT = 


G 


ex{0) 


-1 

0 

^ph(0) 

0 

(0) - 

0 

0 

^-1 

'^ex(0) 


0 

0 

0 


t') and 



0 

0 

0 

^-1 

^pKo). 


(50) 


j(x, t; x', t') are the 
inverse exciton and photon Green’s functions given by 
(fT51) and (1^ . respectively. 

We can rewrite as G”^ = G)(j^ — S, where S 
is a matrix formed by self-energies. Multiplying this 
equation from the left by G(o) and from the right by 
G~^ we obtain the Dyson equation G = G(o)-l-G(o)SG. 
Some Green’s functions from G are presented in Fig. 
1. We will only consider the case of a uniform system 
F“(x) = 0. The exciton and photon single-particle 
energies have the form 




^ph 


= eP^ = 


2mex ^ ‘2mph 

Substituting p, Eqs. (fT^ and (1^ into gives 


(51) 


- ?iG"^(p,fc) = 

'E^^- iTwok m/2 Voxl 0 

hn/2 EP^-ihujk 0 0 

FoXo 0 m /2 

0 0 nn/2 EP>^+inujk 

(52) 





+ - 




□□OddddO = -->0=;;:::=idddO 

b 




"{> 


.{> + 




= ->^=0=^ 

d 

FIG. 1. Diagrammatic representation of some Dyson 
equations from G = G(o) -I- G(o)SG. Solid and dotted lines 
represent free-held exciton and photon Green’s functions, 
respectively. Double wiggly line represents exciton-exciton 
interaction counted in the ladder approximation. Circle 
represents the Rabi splitting, (a) Exact exciton Green’s 
function, (b) Exact photon Green’s function, (c) Exact 
anomalous exciton-photon Green’s function, (d) Exact 
anomalous photon-exciton Green’s function. 
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where we introduced the energies as 

^ ^ ^ 2yono = ^ 

-\J{Eo + Vono)^ + n^n^^ + (t + ^°^°) ’ 

_ ^ = + i 

" (t +^°”°) ■ 

Hereafter, we denote by Uq and replace xq by 


B. Spectrum 


The determinant of the matrix G~^(p,w) is zero at 
the poles uj of the matrix G(p,w). Equating the deter¬ 
minant to zero yields 

B'c.* - + Ef + 

+ - VS<EiA = 0- (55) 

Solving Eq. (1551) for we obtain the polariton energy 
spectrum modified by the condensate: 


r^^(EU)2 ^ 1 J^^ex2 ^ ^ 

±i + VSnif 

+ h^^f{^{El^ + Eff-VSnt)\. (56) 


In the absence of the condensate no = 0, it follows from 
Eqs. (IMl), (El that E^^ = -b Eq - /r, = ef - fi, 

and equation (l56l) gives 


^ ((^0 + + sf + ^ + 

-2/r(Eo+e“+£P'‘)) 

X y(Eo+£--eP'‘)" + ?i"H2, (57) 

which is equivalent to 

;i4^’^) = i(Eo+£“+£^'‘) 

± - p=4^’^) - p. 

(58) 


One can see that the spectrum turns into the energy 
dispersion of the lower and upper polaritons. 

In the absence of the Rabi splitting we have 


i + Ef - 

±l{-E;^^ + EP>^" + Vo^nl), (59) 

which gives the photon energy spectrum 

hujp = ±EP^ = ±{eP^ - m) (60) 

and the exciton energy spectrum modified by the con¬ 
densate 


hu;p = - V^nl = = ±e^. 

(61) 

The equation El is the Bogoliubov dispersion law. 
Here we assume p = Eq -b Von-o as discussed in Section 

4. 

For small momenta p —> 0 the lower branch takes the 
phonon-like form 

= VsP, (62) 

where Vs is the sound velocity. Assuming that a = 
fELph/fELex = 0, VoTio <C hn^ which is commonly encoun¬ 
tered in practice, the velocity can be estimated as fol¬ 
lows: 


/ Vbno 

E^ + imph 

At p —>■ 0 the upper branch has the form 


(63) 


= A -b 


P 


2m, 


(U)-' 


where 


(64) 


A = \ E^ + h^n-^E- 


1-b 


3Eo 


^E^ + n^n\ 


Vono (65) 


in the limit Vono <C hfl, and 



in the approximation a = mph/m^x = 0, Vono <C hn. 

In the limit of large momenta p —>■ (X) the lower polari¬ 
ton branch turns into the dispersion law of free excitons 


^ 

P 2me 


1 

+ 2 


{Eq + Eq^o)^ + -b Eg -b 3Vono^ . (67) 





























If hn = 0, the square root should be taken with the 
negative sign; this corresponds to ^ = Eq + Vbno- 
The upper branch for p ^ oo takes the similar form 


G-i = WG-i is a non-Hermitian operator with re¬ 
spect to the scalar product (a, 6) = Let us 

define the new scalar product as 




+ — (L'O + Lono)^ + — (Eq + Vono)^ . (68) 


The outlined properties of the polariton spectrum mod¬ 
ified by the condensate are shown in Fig. 2. We see 
that the energy and momentum of the polariton reso¬ 
nance increase with the condensate density. This resem¬ 
bles a positive detuning between the exciton and photon 
modes. 


C. Polariton basis 

Let us now diagonalize the matrix G~^. The diag- 
onalization does not reduce to straightforward deter¬ 
mination of the eigenvalues because of opposite signs 
in front of ihujk in the diagonal elements of the ma¬ 
trix (15^ . An option is to consider the generalized 
eigenvalue problem G“^^ = (fiio — ihuJk) W^, where 
W = diag{l, 1, —1, —1) (see, for example^). However, 
we suggest a different approach: to multiply the two 
lower rows of G”^ by — f and thus obtain the stan¬ 
dard eigenvalue problem G“^^ = (huj — ihujk) C, where 



0 1.0 2.0 
P 


FIG. 2. (Color online) Polariton energy spectrum modified 
by the condensate. The momentum is expressed in the 
units of ^^niphTiQ,, the energy is expressed in the units of 
TiQ,. (a)Vbno = 0, (b)Fono = O.llifl, (c)Foa.o = 0.25?in, 
(d)Vbno = QATiQ., {e)Vono = O.Shfi. Inset: phonon-like 
spectra in the vicinity of zero momentum. 


((a,6)) = (a,W&). (69) 


The Hermitian conjugation with respect to the old and 
new scalar products is given by the formulae 

(Oa, b) = (a, O+b), ((Oa, b)) = {{a, 0®&)). (70) 


One can obtain the following relationship between the 
operators O® and 0+: 

o® = WO+W. (71) 

The new scalar product is chosen in such a way that the 
equation (G~^)® = G“^ holds. Therefore the operator 
G“^ is Hermitian with respect to the scalar product (1691) 
and its eigenvalues are real: G“^^ = A^, A = A , if the 
norm of the vector ^ is nonzero: ((C) 0) 7^ 0. 

Let us perform the transformation from the exciton- 
photon basis (x,'0)X)V’) to the polariton basis 
(P(^),P(^),P(^),P('^)), where P(^) and P(^) describe 
the lower and upper polaritons, respectively. We shall 
carry out the canonical transformation 


p{L,U) ^ p,{L,U)y 

p 'Xp Ap 




-XL.U)^ 




'xp 'ap + V'p 

(72) 

where unknown 

coefficients to be found. 

It is clear that 


p{L,u) ^ pL,u)- 

P 'Xp Ap 




4t’"^Ap + 4r^p 


(73) 

and therefore the two bases are related to each other by 
the equation 


pi^) 


r w 

(L) 

^(.L) 

PXp 

-{L)l 

U , 

'ipp 

p{U) 


(U) 

dxj 

{U) 

Pxp 

-{U) 
u, 

'//’p 

p(-^) 

^p 


Pxp 

'^p 

-(L) 

dxj 

-(L) 

p(^) 

^p . 


L^xp 

{U) 

u, 

'Ipp 

AU) 

/Xp 

-(U) 



Xp 


tpp 


Xp 


-V'p- 


(74) 


Let us denote the matrix in the right-hand side by B. 

The diagonalization of a Hermitian matrix implies the 
use of a unitary matrix, which columns are normalized 
eigenvectors of the Hermitian matrix. Let us assume 
that the columns of the matrix are the normalized 
eigenvectors of the operator G~^, i.e. we have 

G-iB+=B+Gd\ (75) 


where 



diag 




ifLUJk, — ihoJk, 

- ihuJk, — ihujk^ . (76) 
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The columns of the matrix B+ are normalized with re¬ 
spect to the scalar product ((•)), i.e. their components 
satisfy the equation 




(L,U) 


y{L,U) 

Xp 


iL,U) 


= 1. (77) 


The components of the eigenvectors can be found from 
the equations = {Tilo — iTitOk) C- We obtain the 

following relationships between the coefficients of the 
transformation: 


iL,U) 


hn/2 






,XL,U} ^ 
Xp 






X Vonoe 


-2i<j, (L,U) 

'Xp > 


(L,U) 


hVL/2 


X Vonoe-^^^ (78) 


Using ([77]) and (l78|) we can find 
and 4p’^^- 

We present another method for the calculation of the 
transformation’s coefficients. From Eq. (EH) we obtain 


G"^ = B+Gd\B+)~i = B+Gd^(B+)® 

= B+Gj^WBW. (79) 

If the matrix G“^ = WG“’^ corresponded to the diag¬ 
onal matrix Gj”^ = WG^^, we would write 

G"^ = B'^Gd’^WBW. (80) 


through the components of B. Thus, for example, 


Gii(p, fc) = 




2 


2 


Vxp 

-7TT + - 

'^Xp 



—iTiujk + fiLUp ihujk + ^ojp 




2 


2 


%p 

L 

^Xp 



—ihujk + 


ifiuik + fiuj] 


(U)-' 


G22(P,*) = 



JL) 

2 


2 



-7TT + - 

‘IPp 



(83) 




2 


2 



L 

Ipp 



—ifiujk + ihuJk + 


(84) 


On the other hand, we can find the matrix G directly 
by inverting the matrix G~^ from Eq. (I52F Thus, for 
Gii and G 22 we have 


Gii(p,fc) = - 


D L 




(EPf 


t.2 2 

h uJk 


ihiOk + El^EP^^-EP^^ ^ 


(85) 


G22(p,A;) = 


h 

~ D 

EPj^n^tol + (e -2 + 

(hn) 

+ 4 

- - ihujk + e;-^ep>^ 


pexi^^) pphy 2 2 

4 *^ono 


( 86 ) 


From this equation, it follows that 


where 


G = B+GdB, 


where 


(81) 


D = {ihuJk — hu}p){ihu}i~ — fiijJp) 

X {ihuJk + fiujp){ihoJk + ). (87) 


Gd = diag 


^ —ihuik + —ihuJk + 


ihuJk -f huj^^ ihuJk + hu!^'^ 


(82) 


Multiplying the matrices Gd, B and B+ according 
to we express the components of the matrix G 


Bringing these equations to the form similar to that 
of Eqs. (1551) and (1551) . we obtain the coefficients 

(L,U) {L,U) , (L,U) 

^V-p ’ Xp : and . 

Both methods lead to the following expressions for 
the coefficients of the transformation from the basis of 


excitons and photons to the basis of lower and upper 
polaritons in the presence of the condensate: 
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p{L,U} ^ gi? 


(L,U) 




\ 




huj. 


(,L,U) _ j^ph 


- Ef'^ Xv + 

+ Ef'^ Xp - 


+Vono 


h.(jj 


i^iU) _ 


\ h^nyA - (ei’^ + nJp'^^^ 


■ ( 88 ) 


Here are the phases which up to an arbitrary 

phase ^ satisfy the equation — -k/2 = ^. 

In a sense, these transformations incorporate both the 
Hopfield transformation for polaritons and the Bogoli- 
ubov transformation for weakly-interacting Bose gas. 

For hEl = 0 the expression (1551) in the case of lower 
polariton yields 


p(.L) 

P 


e-p- + Vono + e^ 

V— 

+ ' V 2cg 


(89) 


The coefficients of Xp n Xp are the coefficients Up and Wp 
of the Bogoliubov transformation. An upper polariton 
turns into a photon: = E’^ijjp. For Vbno = 0 the 

expression (1551) yields 


p(.L,U) 

p 



Here the phase ^ is assumed to be zero. The coefficients 
of Xp and '(/'p are the Hopfield coefficients. 

In Fig. 3 we plot the momentum dependences of the 
calculated coefficients. It is seen that the coefficients for 
the lower polariton similar to the Bogoliubov transfor¬ 
mation coefficients have a singularity at zero momen¬ 
tum. This is due to the fact that, by virtue of Eq. (EH), 
huj^^ in the denominator of Eq. (1551) behaves like Vsp. 
In addition, the presence of the condensate shifts the 
balance between the excitons and photons to the pho¬ 
tons in both polariton branches. Apparently, it is a con¬ 
sequence of the positive detuning between the exciton 
and photon modes discussed in Section 5.2. It is worth 
pointing out that exciton-photon interaction causes an 
extremely slow decay of the coefficient . 


D. Occupation numbers 


We can find the noncondensate exciton density taking 
into account that = —Gii(x,r,x', t+). Using Eq. 



FIG. 3. (Color online) The transformation coefficients 
from excitons and photons to lower and upper polaritons as 
functions of momentum. The momentum is expressed in 
the units of yniphTiO. (a)Vbuo = 0, (b)Vbuo = O.lTiQ, 
(c)Vbno = 0.25?ifi, (d)Vono = O.dhfi, (e)Vbno = 0.5hfi. 
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we obtain 


Similarly, for photons one has 


’^'= - E G'ii(p,A) = 


= -E 




2 


2 


%p 

-7TT -f - 

^Xp 



p/o.fe I ifiuJk + TiWp^^ 




2 


%p 



^Xp 


—ihujk + huj^^ ihujk + huj^'^ 


y (91) 


To evaluate the sum over Matsubara frequencies, we use 
the formula 


1 

lim — > -;- —— 

Tifi ^ iojk — (e — 

and find that 


1 


j/3(e-/i) _ 2^ 


(92) 



(93) 


Substituting the explicit expressions for the coefficients 
^(L.C/) (L,U) (L,U) (L,U) 

’ '^Xp > ^V>p 

as follows: 


^Xp ' I ''I'jb b ^Xp b b we can write this equation 




V 


(94) 


V p#o 


p#o 


where the first term contains the sum over exciton oc¬ 
cupation numbers of the lower polariton branch 


= - 


-EP'^^ + n^uj^ - (nil)2/4 


2 

2, ,l2 


) 


-h 


-E^^EP^ + E^^h^u}^ + EP'^{hn)^/4: 


f-2 r2 ^2 r/2 

n uj^ -n uj^ 

1 /I 


hojp V 2 


‘'p -1 


(95) 


and the second sum is over exciton occupation numbers 
of the upper polariton branch 


^exC/ 

«p 


, -EP’^'^ + - {nnf/A 

2{n^u:^^-n^u2f) 

+ EP'^ihny/4: 


f -2 t 2 f .2 jt 2 

n u)„ —hwii 


TiUjU \2 


"p -1 


■ (96) 


^ph 


' = -G22(x,T,x',T+) = - ^ G22(p,fc) = 


p^0,/c 


= -E 




2 


2 



-rw -f - 

Ipp 



p/o,fc —iTEJk + ^iWp ihuJk + 




2 







2 


i/’p 



—ifiuJk E hujy^ ifiuJk + huj^'^ 


= 



(97) 


which leads to 


-^ph 


' = V IE "p“' +(»«) 


i,P#0 


p#o 


where 


^ -E;^^ + h^u;y - ihnr/A + V,^nl 
2 

-El^'^EP^ + EP^h^uj^ + E^^{hn)'^/4: + EP'^V^n, 


2 t2 ^2 r/2 

^p 

1 /I 1 


huj^ \2 


(99) 


^pnu' ^ -E;^^ + n^ojf-{nn)y4 + VoVo 

^ex2 -t EP’^h^ujy + E^^{nnf/4: + EP’^V^nl 


f -2 t 2 ^2 r /2 

n uj^ -n uj^ 


1 




"p -1 


( 100 ) 


Here and are the photon occupation num¬ 

bers of the lower and upper polariton branches. 

Let us now analyze the results. In the absence of the 
Rabi splitting the first term in Eq. ((551) equals 1/2 and 
the brackets’ coefficient in the second term turns into 
(e“ -I- Vono)/{hujp), where hojp is given by Eq. (I5T]) . 
Both terms in Eq. (1551) are equal to zero, i.e. there 
are no excitons in the upper branch. As a result, for 
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excitons we obtain the expression 


nf = = 

gp + Vonp — hujp £p + Vonp 1 I'mi'i 

2?iWp hojp — 1 

which is the standard momentum distribution of non¬ 
condensate particles in the Bogoliubov model. 

Both terms in (IMl) are equal to zero, i.e. there are no 
photons in the lower polaritoii branch. The first term 
in (llOOl) equals —1/2 and the brackets’ coefficient equals 
1. Thus for photons we obtain 


nf' = nf -k nf - 

P P ' P _ 2^ 


( 102 ) 


further behavior of the exciton iioiicondensate distribu¬ 
tion is shown in the inset of Fig. 4,a. It is seen that 
the distribution drops off at very large momenta. It is 
worth noting that the population of the upper polari- 
ton branch is high (comparable with that of the low 
polariton branch). This population increases with the 
condensate density. The maximum of the population 
corresponds to the region of the polariton resonance and 
shifts to large momenta with the growth of the conden¬ 
sate density. 


VI. COUPLED 

BOGOLIUBOV-DE-GENNES-LIKE EQUATIONS 


which is the standard distribution of particles in the 
ideal Bose gas. In Fig. 4 we plot the noncondensate 
exciton and photon distributions at zero temperature. 
At small momenta the exciton distribution is linearly 
decreasing starting from some finite value: 


In order to investigate the inhomogeneous case we 
introduce the operators 

= + (104) 


27rpn 


exL' 

P 


TT 

2 


TT Vonp A _ 

4 Us V ^ / 



VoUq 


P, 


(103) 


^ph — 




2m. 


ph 


and solve the following eigenvalue problem: 


(105) 


but then it starts increasing with growth of momen¬ 
tum. Such behavior is a consequence of very slow drop 
of the coefficient (see the previous subsection). The 






FIG. 4. (Color online) Exciton and photon momentum 
distributions in the lower and upper polariton branches. 
Momentum is expressed in the units of 
(a)Vbuo = 0, (b)Vbuo = O.lhfl, (c)Vbuo = 0.25hfl, 
(d)Vbno = QATiQ., (e)Vouo = 0.5?ifl. The inset shows the 
exciton momentum distribution in the lower polariton 
branch at large momenta (of the order of 
^Jmexlrriph ~ 10 ^). 


'iCeo; + 2Uo|xo(x)P hQ,/2 Uo(xo(x))^ 0 

nVtl2 Kph 0 0 

Vb(xo(x))^ 0 Aea: + 2Vb|xo(x)P hn/2 

0 0 nn /2 kph. 


■<"(x)' 


0 

0 

0 

1—1 


■<"(x)' 

uf(x) 


0 10 0 


<{^) 

<"(x) 

— 

0 0-10 


u-(x) 

uf(x) 


0 0 0 -1. 


uf(x) 


(106) 


Thus, we obtain a set of equations which is an analog of 
the Bogoliubov - de Gennes system of equations and de¬ 
scribes excitations in the system of coupled condensates 
of photons and excitons: 


2„2 


2me 


+ U“(x) - p + 2Uono(x) <"(x) 


hn 


-hVb(xo(x))^u“(x) -h —<''(x) = hujnul^ix), 
<'‘(x) + ^Mr(x) = ;iw„<'‘(x). 


2^2 


V 

2me 


-h U“(x) - p -h 2Vono{:>c) ) u“(x) 

nn 


+Uo(xo(x))^ur(x) + —uf (x) = -a.;„ur(x), 

(■£$ - (x) + ^C(x) = (x). 

(107) 
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Similarly to Eq. (1931) , one can write an expression for 
noncondensate excitons 


VII. CONCLUSIONS 


n/0 


and for photons 


(x)|VK-(x)|^ 

e/3‘-n _ 1 


+ k; 



(108) 



(109) 

Functions u‘^{x), ^“(x), u^{-x) and w^^(x) are nor¬ 


malized as 


J c^x(|u“(x)|V|<'‘(x)|^- 

|ur(x)|"-|<^x)f) =1. (110) 


In conclusion, the central result of the work is the 
canonical transformation from the basis of excitons and 
photons to the basis of lower and upper polaritons in the 
presence of the condensate. The transformation coeffi¬ 
cients we obtained can be called the Hopfield coefficients 
renormalized by the condensate. 

The path integral method we used has essentially 
more abilities than were demonstrated in the work. The 
method opens an effective way to analyze the formation 
and dynamics of a non-equilibrium condensate. The 
two-component approach based on the path integral 
method leads to the two-component diagram technique, 
which was briefly presented (Fig. 1), but not applied. 
All this aspects will be discussed in future works. 


This equation is analogous to Eq. (EZD and can be ob¬ 
tained in the same manner. 

The solutions of the set of equations (11071) are out of 
the scope of the present discussion and will be published 
elsewhere. 
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